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Abstract 

A method is presented to bend a thin massive hne when the curva- 
ture is small The procedure is apphed to a homogeneous thin bar with 
two types of curvatures. One of them mimics a galactic bar with two 
spiral arms at its tips. It is showed that if the bending function is a 
linear combination of Legendre polynomials, then the bent potential is 
an exact solution of the Laplace equation. A transformation is applied 
on the thin bent bars to generate three-dimensional potential-density 
pairs without singularities. The potentials of the thin bent bars are 
also used to generate non-axisymmetric planar distributions of matter. 
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1 Introduction 

Bars are a common self-gravitating structure present in disc galaxies. About 
50 per cent of such galaxies are stongly or weakly barred, including our Milky 
Way [U [2] ; see also the classification of galaxies by [3] and the fraction of 
barred galaxies discussed by [HEIE]. Galactic bars are triaxial systems, and 
constructing analytical triaxial potential-density pairs is a difficult task. The 
only exact, self-consistent models of bars were constructed by Freeman [7], 
but they have some unrealistic features for barred systems. As alternatives, 
galactic bars have been approximately modelled as homogeneous ellipsoids 
[U [9] or inhomogeneous prolate spheroids \10\ [TTl [T2| I13j . In these works. 
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the inhomogeneous bar has been represented by a Ferrers elhpsoid |14] . 
which has a finite length and represents many features of galactic bars rather 
well. Long &: Murali [TS] found simple analytical potential density-pairs for 
prolate and triaxial bars that can all be expressed in terms of elementary 
functions. One of their model of bar was used in hydrodynamic simulations 
[TBI [T71 [TEl dHl [20] to understand the response of a gaseous disc to the 
imposition of non-axisymmetric bar potentials. 

Until recently, an unclear issue was the connection between bars and 
grand-design spirals. Many barred galaxies have spiral arms that appear to 
emerge from the tips of the bar (see, e.g., NGC 1300; [21j, p. 525). How- 
ever, [22] present evidences that the pattern speed of the spirals are much 
lower than the pattern speed of the bar, so the spiral cannot be driven 
directly by the bar. On the other hand, there are other observational evi- 
dences that bars and spiral arms are correlated [231 IMl HSj [26]. Salo et al. 
[27) recently investigated the relation between bar forcing and spiral den- 
sity amplitudes for over 100 barred galaxies, and found that there exists a 
significant statistical correlation. Furthermore, hydrodynamics simulations 
of the response of a gaseous disc to the imposition of a non-axisymmetric 
bar [28l[Ml[3Ql[3ll[IZl[32l[33l[Sl[35l[l9] have shown that the symmetric 
two-armed spirals in barred galaxies are driven by the gravitational torques 
of the bar. 

The system bar-l-spiral arms may be viewed as a bar with bended ends, 
although we are not aware that such an interpretation has been proposed 
so far. It would be interesting to have simple, analytical models for such a 
gravitating system. In this work we propose a method to obtain potential- 
density pairs for thin and for 'softened' bent bars. In Section [21 we present 
a procedure to bend a thin massive line and to calculate its potential. The 
idea is to consider a slight curvature and expand the potential with respect 
to a small parameter. In Section [3| this formalism will be particularized to a 
thin bar with constant linear density. It will be shown that if the 'bending 
function' can be written in terms of Legendre polynomials, then the poten- 
tial of the bent bar will be an exact solution of the Laplace equation. Two 
examples of deformed bars will be discussed. By using a suitable transfor- 
mation, the thin bent bars are then 'softened' to generate three-dimensional 
potential-density pairs without singularities. In Section [H we present non- 
axisymmetric potential density-pairs that represent planar distributions of 
matter constructed from the two potentials of bent bars discussed in Section 
[3l These planar potential-density pairs are found by using a method first 
proposed by Kuzmin [36j . The discussion of the results is left to Section \5\ 
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2 Bent massive lines 



In this section we present a procedure to bend a thin massive hne. The 
curvature is supposed to be small, and the potential of the bent system 
will be obtained from an expansion of a small parameter. The undeformed 
massive line will be described by the parametric equations x' = xoit), y' = 
yo(i), z' = zo(0) with t G [ii,t2]- The gravitational potential $(r) of this 
line can be expressed as 

where G is the gravitational constant, s is the arc length and A(s) is the 
linear density. In terms of the parameter t, equation ([1]) is rewritten as 

* = -a r '-P^, (2) 

Jti \r-r'\ 

where the dot represents derivative with respect to t. 
We shall represent a slight curvature of the line as 

r' = ro(t) + eri(t), (3) 

where ro(t) is the position of the undeformed line, e is a small dimensionless 
parameter and the particular form of the curvature will be determined by 
the function ri(t). We have, up to first order in e, 

|r'| = |ro| +£— — , (4) 
|ro| 

1 1 , ^ ri • (r-rp) 

I 7T = I f + ^n 13~- ^ 

|r — r| |r — ro| |r — ro|'^ 

Thus, up to first order in e, the potential ([2]) can be written as 
^ ^ _^ j-'^ X{t)\iro\dt Z-*^ A(t)ro-ridt f'^ A(t)|ro|ri • (r - tq) dt 



rol Jti |ro||r-ro| Jt^ |r - roP 

(6) 

3 Bent bars 

Now we consider the particular case of a bent bar with constant linear den- 
sity. The undeformed bar located symmetrically on the z-axis with length 2a 
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and linear density Aq will be parameterized by tq = (0, 0, t), with t € \—a, a 
For simplicity, the 'bending function' ri{t) will be chosen as ri 
With these assumptions, the potential ([6]) reduces to 



$ = $0 + e^i, 



o,/(t),o). 

(7) 



where 



$1 



-GAo 
-GAo 



dt 



2-2 _j_ y2 _j_ 



yf{t)dt 



x"^ + y"^ + {z — t) 



3/2 ■ 



At this point it is convenient to relate the potentials 
found by Letelier 



Qn{u)Pn{v) 



Pn{t/a)dt 
R^ + {z- tf 



(8) 

to an identity 
(9) 



where P„ and Qn are, respectively, the Legendre polynomials and the Leg- 
endre functions of the second kind, and (u, v) are the spheroidal coordinates 
related to the cylindrical coordinates (i?, z) through 



{Ri+R2)/{2a), 



Ri = JR'^ + {z + ay 



v = {Ri-R2)/{2a), 
R2 = \/r'^ + {z ' 



a 



(10) 
(11) 



with u > 1 and —l<v<l. The physical interpretation of ([9]) is in terms 
of a bar with linear density proportional to a Legendre polynomial, whose 
potential corresponds to a multipole term that arises in a multipolar solution 
of the Einstein equations. The potential $0 in dS]) can thus be identified with 



$0 = -2GXoQo{u)Po{v) = -2GXoQo{u). 
Furthermore, we note that <I>i can be rewritten as 

yf{t)dt d 



(12) 



-GAn 



+ y2 + (z — t)^ 



3/2 



GAn^ 



f{t)dt 



dy J -a . 1^2 + y2 + (2 _ if 

(13) 

If the function f{t) can be expressed as a linear combination of Legendre 
polynomials, then the last integral in p3|) may be directly integrated via the 
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Figure 1: The shape of the bent bar for (a) the function f{t) = t'^/a, and 
(b) the function f{t) = t^jc?. 

identity ([9]). This also shows that the potential ^I^i will be then a solution 
of the Laplace equation, thus the potential ([7]) of the bent bar will also be 
an exact solution of the Laplace equation. 

As examples, we will calculate the potential-density pairs for bars with 
two choices for the 'bending function' ri = (0,/(t),0). By (l3|), the para- 
metric equations of the deformed bar are x' = 0, y' = ef{t) and z' = t, 
t € [—a, a], which may be rewritten as y' = ef{z'), —a < z' < a. We shall 
take as first example f{t) = /a and then f{t) = /a?' . The shape of the 
bent bar in each case is schematically depicted in Figs[l]^a) and (b), respec- 
tively. The first shape does not represent necessarily a galactic bar, but we 
chose this bending function for its simplicity. The second shape mimics a 
bar with two spiral arms. 

In terms of Legendre polynomials, we have 

t^^pJl)Ajl), (14) 

a 3 \a J 3 \a J 

5 \a J 5 \a J 
It is then straightforward to calculate the potential (jl3p using the identity 
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([9]), and we state the final result 



(a) 



GAo In 



+ e 



+- 



CApj/ / z-a + R2 
a \z + a + Ri 

{Ri + R2) + z (x^ + y^ + z^) {Ri - R2) 



z — a + R2 
z + a + Ri 



GAo In 



(x2 + y2) R^R^ 

z — a + R2' 



(16) 



X <; 3z (x^ + y'^) RiR2ln 



z + a + Ri 
z — a + R2 
z + a + Ri 



GXoy 



+ e 



+ az (5a;^ + 5y^ 



a? (x2 + y2) _R^_R2 



- [a^ (x^ + y') + (x^ + y2 + z^) {2x^ + 2y' - z^)] } , (17) 



where i?i = \J x"^ + y"^ + {z + a) , R2 = \J x"^ + y"^ + {z — a) and the sub- 
scripts (a) and (6) refer to the potentials calculated with the functions 
f{t) = t^/a and f{t) = /a^, respectively. The potential (fTHl) remains 
invariant under the transformations x — )■ —x or z — >■ — z, whereas the 
potential (llTh remains invariant under the transformations x — >■ — x or 

y ^ -y,z ^ -z. 



3.1 Softened bent bars 



The potentials (I16p - (|17p become singular along the bent thin bars. In order 
to get physically more realistic potentials it is convenient to 'soften' them. 
A simple way to achieve this is by means of a Plummer-like transformation 
|15j . In our examples we apply a transformation x^ + — : 



x^ + y^ + h\ 

where 6 > is a 'softening' parameter. The corresponding three-dimensional 
mass density distribution is obtained from the Poisson equation in Carte- 
sian coordinates, V^^>(x,y,z) = AirGp. From the potential (jl6p . after the 
transformation, we obtain the pair 



(a) 



GAo In 



z — a + 7^2 
z + a + TZi 



+ e 



GXoy 



In 



z — a + 7^2 
z + a + TZi 



+ 



a (x^ + j/^ + fe^ - z^) (7^l + 7^2) + z (x^ + + fe^ + z^) {Ui - 712) 

(x2 + 2/2 + 52)7^^7^2 



Pia) = P0 + ePl, 



(18) 
(19) 



where 



Po 



-n\ {z - a) 



{7^i {z + a) [3 {x"^ + y V 6^) + 2 (z + a)^ 



3 (x^ + y^ + b^)+2{z- a)2] } , (20) 

- |6z (x^ + y2 + 62^3 _ 
47ra(x2+y2 + 52)3 ^5^5 \ V ; I 2 1; 

-Tef (z - a) Isz'^ {z -a)^ + 2{z- af (x^ + + 6^) (Hz^ - 2az + a^) 
+5 (x^ + y^ + 6^)^ (4z2 - 2az + a^)] + (z + a) [sz^ (z + a)^ 
+2 {z + a)2 (x^ + y2 + 6^) (11^2 ^ + a^) 

+5 (x^ + y2 + 62)^ (4z2 + 2az + a^)] } , (21) 



Pi 



+ y2 + 62 + (z + a) , and 7^2 = V + y^ + 62 + 



where TZi = 

Note that the first term of the potential (jlSp is the same as the potential of 
the prolate bar of |T5] with the replacements z ^ x and x^ + y^ ^ y'^ + = 
R\ 

Figs[2]^a)-(c) display some curves of constant density of the dimensionless 
mass density p(^a) = P(a)/{^o/(i'^), equation (fT9]l . in the three orthogonal 
planes, with parameters b/a = 0.25 and e = 0.15. The bending of the bar 
is apparent in the y — z plane in Fig. [2|^a). The density contours near the 
origin are bent in the same manner as the thin line of Fig. [TJ^a). There are 
two points of maximum that are shifted to the right of the z-axis, and two 
points of minimum appear on the opposite side. At these points the mass 
density may become negative as e is increased; we found that for b/a = 0.25 
the density is non- negative for e ^ 0.16. As the parameter b/a is increased 
for a fixed value of e, the density contours become less elongated, the two 
points of minimum may disappear, and the two points of maximum move 
towards the origin. The contours in the x — z plane preserve axial symmetry. 
Some level curves of the potential = <I'(a)/(GAo), equation (fT8]l . with 
parameters b/a = 0.25 and e = 0.15, are displayed in Figs [3)^a)-(c) . In the 
y — z plane the potential has only one minimum on z = and y = ym , which 
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is the zero of 



_^ £ / V^m + h'^ + a? 



{Vrn + b"^) Vvl + b'^ + a^ « \y^yl + b^ + a'^ + 



2e 



{vL + fe^) {via + + a' 



2^3/2 



0. (22) 



For the above-mentioned values of parameters the minimum occurs at ymlo' ~ 
0.011. The components of the force '^(a) = —'^^(a), corresponding to the 
potential p^ . are listed in Appendix Rl 

Applying the same transformation on the potential (I17p results in the 
potential-density pair: 



GAn In 



z — a + IZ2 
z + a + IZi 



a? (x2 + y2 + 52)7^^7^2 



+az (5x2 ^ 5^2 ^ 5f,2 _ ^2^ (7^^ ^ _ _ ^^2 ^^2 ^ y2 ^ j,2^ 



where 



Pi 



+ {x^ + y^ + b^ + z^) (2x2 + 22/2 + 26^ _ z^)] } 



(23) 
(24) 



iz^ (z + af 



W (x2 +y2 + ^2)3 7^57^5 

+2z (z + a)^ (x2 + y2 + ^^3^2 _^ ^ 3^2^ 
+15z (z + a) (x2 + y2 + ^2^2 ^^^2 ^ + a^) 
+ (x2 + y2 + 52^3 ^^^^2 ^ ^Q^^ ^ 5^2^ + 2 (x2 + y2 + 52^4 

-7^f (z - af + 2z{z- af (x^ + + ^2^ (^^^^2 _ _^ 3^2^ 
+15z {z - a) (x2 + y2 + ^2^2 ^^^2 _ + a"^) 

+ (^2 ^y2^ ^2)3 ^^^^2 _ 10^^ ^ 5^2) + 2 (x2 + y2 + 62)4J | ^ ^^5) 

and po is given by ()20p . Some level curves in the three orthogonal planes of 
the mass density /3(b) = p(j^^/ a?), equation (p^ . are shown in FigsHl^a)- 
(c) with parameters b/a = 0.25 and e = 0.15. The density contours near 
the origin are bent in the same manner as the thin line of Fig. [^b). Two 
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points of maximum are present near the z-axis and two points of minimum 
are located antisymmetrically. The mass density may also become negative 
at these points with increasing values of e, which can rise up to e ~ 0.18 
for h/a = 0.25. When the parameter h/a\s increased with fixed e, the point 
of minimum may disappear, the points of maximum approach the origin 
and the isodensity contours become more circular, in the same manner as in 
the previous case. Figs [5^a)-(c) display isopotential curves of the potential 
= $(6)/(GAo), equation ([23j) . with parameters h/a = 0.25 and e = 0.15. 
Here the minimum of potential is located at the origin. The components of 
the force F(f,) = — V<I>(;,) corresponding to the potential ([23]) . are also given 
in Appendix Rl 



4 Non-axisymmetric thin distributions of matter 

Non-axisymmetric thin distributions of matter are of interest since many 
galaxies exhibit asymmetries in their discs (see e.g. [3_8J [39] for observational 
evidences of asymmetries in disc galaxies). The potential of a thin bent bar 
may be used to generate non-axisymmetric thin distributions of matter. The 
procedure is similar to the one proposed by Kuzmin [36] to find the gravity 
field of an axisymmetric disc. A source of gravitational field is placed below 
a plane 2 = 0. Above the plane this gives a solution of the Laplace equation. 
The solution valid for z > is then reflected with respect to z = so as 
to give a symmetrical solution of the Laplace equation above and below the 
plane. This introduces a discontinuity in the normal derivative on z = 0, 
which by the Poisson equation results in a surface density. Kuzmin applied 
his procedure on a point mass. In general relativity, where the Schwarzschild 
solution in Weyl's metric is represented by a finite bar with constant density, 
Kuzmin's method has been used to construct 'generalized Schwarzschild' 
discs [:40j. Some other examples of general relativistic discs generated by 
Kuzmin's method can be found in |4H H2l H3l I44j . 

We will place the thin bent bar along the x-axis and then displace it at a 
distance c > below the plane z = 0, and make the reflection with respect 
to that plane. This is equivalent to first make the replacements z — > x and 
x — )• 2 in the potentials p6 p - ()17p . and further apply the transformation 
z — )■ c+ lz|. By Poisson equation, the surface density a on z; = is given by 
[21J 

1 



27rG dz 



(26) 
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For the transformed potentials p^ " (fT7|) we get, respectively, 



_ Aqc [{x + a) R2 - (x - a) Ri] Xocy 
""^"^ " 27r (y2 + c2) R1R2 ^27ro(y2 + c2)2R3R3 

X |r? 2x {x^ + + c^)^ - a (y^ + c^) (a^ - 3ax + Gx^) 

-2ax2 (a^ - 3ax + 3x^)] - R^ 2x {x^ + + c^f 

+a {y^ + c^) (a^ + 3ax + 6x^) + 2ax^ (a^ + Sax + 3x^)] } , (27) 
_ Aqc [(x + a) R2 - (x - a) Ri] Xocy 
""^'^ ~ 27r (y2 + c2)RiR2 ^"2^2 (y2 + c2)2R3R3 

X {r? [2 (x^ + y2 + c2)3 + 3« (y2 + ^2)2 _ 
-3ax (y^ + c^) (a^ - 3ax + 4x^) - 2ax^ (a^ - 3ax + Sx^)] 

-Ri [2 (x^ + y2 ^ ^2^3 ^ 2^ ^ ^2^2 ^ 

+3ax (y^ + c^) (a^ + 3ax + 4x2) ^ 2ax^ (a^ + 3ax + Sx^)] } , (28) 



where Ri = y y2 + c2 + (x + a)^ and R2 = y y2 + c2 + (x — a)^. Note that 
the first term in (I27p and (I28p is the surface density of a flat bar given by 
equation (11) of [15j. Although they obtained the result by convolution of 
a flat Miyamoto-Nagai disc (which is in fact a Kuzmin disc) with a needle 
density, our procedure gives the same result. 

Some contours of the dimensionless surface density a'(a) = cT(a)/(Ao/a) 
are shown in Fig.[6]^a), for c/a = 1 and e = 2, and in Fig.[6]^b) for c/a = 1 and 
e = 0. In both figures the level contours are plot for the same set of values. 
FigsEllc) and (d) display contours of the potential = <I>(a)/(GAo) for the 
same parameters as in Figs[6ja) and (b), respectively. The isopotentials are 
also plot for the same set of values in both figures. Comparing the isodensity 
contours of the flat bent bar in Fig. [B^a) with the flat bar in Fig. M^h), we 
note the asymetry with respect to the x-axis, with two points of minimum 
located symmetrically with respect to the y-axis. For c/a = 1, the surface 
density of the flat bent bar is everywhere non-negative provided e ^ 2.2. 
The isopotential contours of the flat bent bar are shifted along the y-axis. 
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(a) (b) 




-4 -2 2 4 _4 _2 2 4 




-4 -2 2 4 _4 _2 2 4 

x/a x/a 



Figure 6: (a)-(b) Isodensity contours of the surface density a(^a) = 
(7(a) /(Ao/a), equation (i27|l . (a) Parameters: c/a = 1 and e = 2. (b) Param- 
eters: c/a = 1 and e = 0. (c)-(d) Isopotential contours of the transformed 
potential ^(a) = *^(a)/(GAo) for the same parameters as in (a)-(b), respec- 
tively. 
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Figure 7: (a)~(b) Profiles of the surface density = (T(f,)/(Ao/a), equation 
(j27p . along the y-axis, for different values of the parameters, (a) Parameters: 
c/o = 1 and e = 0, 1 and 2. (b) Parameters: e = 2 and c/a = 1, 1.5 and 
2. (c) Profiles of the same surface density along the cc-axis. Parameters: 
c/a = 1, 1.5 and 2. 
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(a) (b) 




-4 -2 2 4 _4 _2 2 4 




-4 -2 2 4 _4 _2 2 4 

x/a x/a 



Figure 8: (a)-(b) Isodensity contours of the surface density = 
(T(f,)/(Ao/a), equation (f28l) . (a) Parameters: c/a = 1 and e = 2. (b) Pa- 
rameters: c/o = 1 and e = 0. (c)-(d) Isopotential contours of the trans- 
formed potential $(f,) = <I>(f,)/(GAo) for the same parameters as in (a)-(b), 
respectively. 
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with a point of minimum at x = and y = ym, which is the zero of 



2aym e ( -s/ym + c^ + a 




(yl + ^yl, + c2 + a2 a \ y^yl + c2 + a2 + a 



yl + c')' + a^ {2yl + 



c 



iVm + C^) iVm + c2 + a 



2^3/2 



0, (29) 



an expression similar to (I22p . For the parameters c/a = 1 and e = 2 we 
have ym/ci ~ 0.376. 

In Figs[7)^a)-(c) we display curves of the surface density ct^^) along the y- 
and x-axes for some other values of the parameters c/a and e. In Fig. El^a) 
the parameter c/a = 1 is kept constant. With increasing values of e, the 
asymmetry of the density profile with respect to y = is enhanced. If we 
maintain e = 2 constant, as in Fig. [Tl^b), the density distribution becomes 
somewhat less concentrated as the values of the parameter c/a increase. Fig. 
nic) shows curves of the surface density along the x-axis for the same values 
of c/a as in Fig. [7l|b). From (|27p . it is seen that for y = the profiles of the 
surface density of the flat bent bar and of the flat undeformed bar are the 
same, since the second term in (j27p vanishes. 

In Figs El^a) and (b) we display some contours of the dimensionless sur- 
face density = a ^p-^ / {Xq / a) , equation ([28]). for parameters c/a = 1 and 
e = 2, and for c/a = 1 and e = 0, respectively. In both figures the level 
contours are plot for the same set of values. Figs[8]|^c) and (d) display con- 
tours of the potential = <I>(b)/(GAo), for the same parameters as in Figs 
[8l|a) and (b), respectively. The isopotentials are also plot for the same set of 
values in both figures. The surface density as well as the potential for this 
model of fiat bent bar are asymmetric with respect to the x- and y-axes, in 
contrast to the potential-density pair of the fiat undeformed bar. For the 
parameter c/a = 1, we find that the surface density of this model fiat bent 
bar is everywhere non-negative provided e ^ 3.3. The potential has a point 
of minimum at the origin. From ()28p it is straightforward to see that the 
surface density profiles along the x- and y-axes are the same for the fiat bent 
and for the fiat undeformed bars. 

In Appendix [B1 we list the components of the force F(„) = — V$(a), 
and of the force F(b) = — V$(b) on the plane z = 0, after applying the 
transformations z — ?> x and x — >■ z, then z — )• c -|- |2:|, on the potentials 

(na-dTTD. 
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5 Discussion 



We presented a method to bend a thin massive hne when the curvature is 
small. The potential of the bent system is obtained from an expansion with 
respect to a small parameter. The procedure was then applied to a homoge- 
neous bar with two examples of 'bending functions'. We showed that if the 
'bending function' can be expressed in terms of Legendre polynomials, then 
the potential of the bent bar is an exact solution of the Laplace equation. 
Potential-density pairs for 'softened' bent bars were constructed by using a 
Plummer-like transformation. The resulting mass density distributions are 
non-negative everywhere for restricted values of the deformation parameter 
£. We also used the potentials of the bent thin bars to construct planar 
distributions of matter without axial symmetry. Furthermore, non-negative 
surface density distributions that are non-symmetric with respect to one or 
both X- and y-axes can be found for restricted values of the deformation 
parameter e. 

We would like to mention that the potential (l6|) of a bent massive line 
contains only terms originated from the expansion with respect to £ up to 
first order (dipole terms). In principle, higher order terms can be included, 
but we found that the explicit expressions, for instance the quadrupole 
terms, are very cumbersome. In particular, the quadrupolar terms can be 
explicitly found using the same algorithm employed for the dipolar terms. 

A future work will be the study of orbits in the potentials of the 'softened' 
bent bars in an axisymmetric background with and without rotation and 
comparing the results with the works with undeformed bars, e.g. |10^ \TT\ 
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A Components of forces for bent bars 



The components of the force F(„) correspondmg to the potential (jlSp . are 
given by 



GXqx [{z - a) 7^l - (z + a) 7^2] 



+ e- 



GXoxy 



|7^? 2z [x^ + y'^ + + z^f -a (x^ + + 6^) (a^ - 3az + Gz^) 



-2az2 (a^ -3az + 3z2)] - 7^: 



2z (x^ + y2 _^ 62 ^ 



+a {x^ + + 6^) (a^ + 3az + Gz^) + 2a2^ (a^ + 3az + 3z^)] } , (30) 

_ GXpy [{z - g) 7^l - + a) 112] GXp / z-a + 7^2 ^ 
^'^^^ (x2 + y2 + 62) 7^l7^2 ^ a {z + a + nj 



GX 







a(x2 + y2 + 62)'7^37^3 

-a (x^ + + 6^ + z^) 



z (x'^ + y'^ + b'^ + z 



2\2 



3^2 (y2 



^2_52)^(^2^y2^^2)2 



-a X 



+ y2 + 6^) [a (x2 + 2y2 + 6^ 



Z X 



+ + 6^)] 



.2^2 (y2 _ ^2 _ f,2) _ 3^) 



+a (x^ + + 6^ + z^) 



/ 2 , 2|i2| 2\2/2 2 

z [x +y +b + z ] [y — x 



3^2 (y2 



2\2 



2 - 6^) + (x2 + y2 + 6' 
(x2 + y2 + 6^) [a (x^ + 2y^ + b^) + z (x^ + 4^^ + 6^)] 



+a^z^ (y2 - x2 - 62) (a + 3z)] } 



(31) 



\a)z 



GXo {112 -n^] 



+ e- 



GXoy 



{ill 



2 (x2+y2 + 62 + z2)2 



7^l7^2 ' ' a (x2 + y2 + 62) 7^f7^3 

+3a (x^ + y^ + 6^) (a + 2z) + 2az (a^ + 3az + 3z^)] 

2 (x^ + y2 + 52 ^ ^2)2 ^ 3^ ^^2 ^ ^2 ^ j,2) _ 3^) 

-2az(a2-3az + 3z2)]}. (32) 
The components of the force F(^) corresponding to the potential (|23p. 
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are given by 



GAox[(- -«)7^l-(^ + «)7^2] 



+ e- 



GX(}.rjj 



{7^? [2 (X^ + y2 + 6^ + ^2)3 ^ 3^ (^2 + ^2 + ^2)2 _ 

-3az (x^ + y'^ + b^) (a^ - 3az + 4^2) - 2az^ (a^ - 2,az + 3^^)] 



2 (x^ + y2 + 52 ^ ^2y3 ^ 3^ / 2 ^ ^2 ^ ^2n2 ^ 



+3az (x^ + + h^) [a^ + Sa^ + 4^2) + 2az^ [a^ + 3a,z + 3z^)] } , (33) 



GAoy [(-z - a) - (z + a) 7^2] GAq-z z - a + 7^2 ^ 



+£- 



(x2 + y2 + ?,2)7^^7^2 



3£- 



z + a + TZi J 

{nl[{x' + y' + b' + z'f 



a2 (x2+y2 + 52)2 ^37^3 
X [2 (x^ + y2 + ft2) (^2 ^ 2y2 + ^2) ^ ^2 ^^2 _ ^2 _ ^2)j 



.2n2 



-az^ (y^ - x^ - 6^) (a2 _ 3^^ ^ 3^2) ^ ^4 (^^2 + ^2 _^ ^: 
-3a (x2 + y2 + ^2)2 ^3^2 + 5^2 + 3^2) _ „ (^2 + 2y2 + ^2)] 
-az (x^ + + 6^) [-3az (4x^ + + 46^) + (x^ + 3^^ + 6^) 



+a2 (7x2 + 102/2 + 762)] 



(x2 + y2 + ^2 ^ ^2 



X [2 (x2 + y2 + ^,2) (^2 ^ 2y2 + ^2) ^ ^2 ^^2 _ ^2 _ ^2)j 

+az3 (y2 - x2 - ^2) (a2 + 3az + 3^2) + (a;2 + y2 + ^,2)2 

+3a (x2 + y2 + 62)2 ^^^^2 ^ 5y2 ^ 3^2) + ^ (^2 ^ 2y2 + ^2)] 
+az (x2 + y2 + 62) [3az (4x2 + 7y2 + 462) + 6^2 (a;2 + 3y2 + ^2) 
+a2 (7x2 + 10y2 + 762)]]|, 



(34) 



GAo (^2-7^l) _3^GAoy^^ 



GXoy 



z — a + TZ2 \ 



'a? (x2 + y2 + 62) 7^37^3 



7^? 



z + a + TZi J 

3Z (x2+y2 + 62 + ^2)2 



-a (X 



+ y2 + 62) (6^2 _^ 3^^ _ 4^2) _^ 3^ ^^2 + y2 _^ ^: 



2\2 



-3az2 (3^2 - 3az + a2)] - 7^i hz {x^ + + + z^ f 

+a (x2 + y2 + ^2) (g^2 _ 3^^ _ 4^2) _ 3^ ^^2 + y2 _^ j^2)2 

+3az2 (3z2^^az + a2)]}. 



(35) 



In all equations, we have TZi = -J x'^ + y'^ + b"^ + {z + a) , and 



^2 = ^Jx^ + y^ + b^ + {z-ay. 

B Components of forces for non-axisymmetric thin 
distributions of matter 

The components of the force F(^) = — V$(a)) the plane z = 0, after 
applying the transformations z — > x and x ^ z, then z — > c + \z\, on the 
potential read 



(a)x 



GXq (R2 — Ri) 



+ e- 



GXoy 



,2\2 



R1R2 ' ~'a(y2 + c2)R3R3 

+3a + c^) (a + 2x) + 2ax (a^ + 3ax + 3x^)] - Rf 2 (x^ + + c^)^ 

+3a + c^) (a - 2x) - 2ax (a^ - Sax + Sx^)] } , (36) 
GXoV [(x - a) Ri - (x + a) R2] GXq / x - a + R2 \ 



+e- 



In 

(y2 + c2)RiR2 a Va^ + a + R-iV 

GXq 



{rI[x (x' + + c')' {y' - c') 



a (y2 + c2)2 r3r3 

-a (x2 + y2 + [3x2 (y2 _ ^2) ^ (y2 ^ ^2)2 

(y^ + c^) [o (2y2 + c^) - X (4y2 + c^)] - a^x^ {y^ - c^) (a - 3x) 

(^2 ^ y2 ^ ^2)2 ^^2 _ ^2) + a (x2 + y2 + c2) [sx^ (y^ _ ^2) + (y2 + ^2)2" 

(y2 + c^) [a (2^^ + c^) + x (4y2 + c^)] + a^x^ (y^ _ (a + 3x)] } . 

(37) 

The components of the force F(t) = — V<I>(fc), on the plane z = 0, after 
applying the transformations z — )• x and x ^ z, then z — )• c + |z|, on the 



-Ri 
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potential (fT7|) . read 

GAo (R2 — Ri) 



{b)x -- 

{r? 



R1R2 



In 



X 



a + R2 \ 



GXoy 



3^ (^2^y2^c2)'-a^-2 



X + a + Ri y a2 (y2 + RfR^ 
(y2 + c^) (6x2 _^ _ 4^2^ _^ 3^ (^y2 _^ ^2^2 
-3ax2 (3x2 _ 3^^ _^ ^2^j _ 3^ ^^2 + y2 _^ ^2^2 
+a (y2 + c^) (6x2 _ 3^^ _ 4^2^ _ 3^ ^^2 + c^)^ + 3ax2 (3x2 j_ 3^^ _^ ^2-^ | ^ 

(38) 



'(% 



GAo?/ [(x - a) Ri - (x + a) R2] „ GAqx 

oe — — m 



(y2+c2)RiR2 



X - a + R2 \ 



X + a 



a? (y2+c2)2R3R3 



{R?[(x2 + y2 + c2)^[2(y2 + c2) (2y2+c2) 



2 1^.2 



+x [y 



ax 



3 /^.2 



y2 - c2) (a2 - 3ax + 3x2) ^ ^4 ^^2 ^ ^2^2 



-ax 



-3a (y2 + c2)^ [x (5y2 + 3c2) - a (2y2 + c2)] 
(y^ + c2) [-3ax (7y2 + 4c2) + 6x2 ^3^2 _^ ^2^ _^ ^2 ^-^^3^2 _^ 7^2^ j 



-Ri 



X 



+ 2/2 + c2)2 [2 (y2 + ^2) (2y2 + ,2^ + ^2 (^2 _ ^2)] 



+ax^ (y2 _ c2) (a2 + 3ax + 3x2) _^ ^4 ^^^2 _^ ^2^2 
+3a (y2 + c2)2 [x (5y2 + 3c2) + a {2y^ + c2)] 
+ax (y2 + c2) [3ax (7y2 + ic^) + 6x2 ^3^2 _^ ^2^ _^ ^2 ^-^Qy2 _^ 7^2^j | _ 



(39) 



In all equations, we have Ri = \/ + + {x + a)"^ and 
R2 = A/y2 + c2 + (x-a)2. 
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